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We show that an impurity embedded in an ionic crystal can give rise to a novel 
kind of local states. These states exist within a polariton gap of a material 
and are a mix of excitations of the crystal, such as phonons or excitons, and 
the transverse electromagnetic field. Electromagnetic component of the states 
along with the corresponding excitations of the material are localized in the 
■ vicinity of an impurity. 
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Impurities or other defects in otherwise ideal crystal lattices have an considerable affect 
on dynamical properties of host crystals. One of the most remarkable dynamical effects 
caused by defects is the occurrence of local states within frequency regions forbidden for 
excitations of the host materials. These states are well known for lattice vibrations |l]-fl; 
excitons magnons |||| and other kinds of excitations in periodic crystals. Such local 

or bounded states had not been found, however, for electromagnetic waves in ideal crystals 
with isolated impurities. The reasons for this are quite clear. For local states to appear 
there must be either bandgaps in the spectrum of excitations of the ideal structure or an 
upper limit of frequencies available for them. Besides that, the excitations and defects must 
be able to interact effectively in order to form a bounded state. Wavelength of light in 
the range of frequencies up to ultraviolet is much greater than the crystal lattice constant 
and characteristics scales of any possible microscopic defects. Therefore, neither of these 
conditions are fulfilled for electromagnetic waves at optical frequencies propagating through 
regular crystals. Bounded photon states were eventually created in so called photonic crys- 
tals flnj. These are periodic arrays of macroscopic pieces of materials with a lattice constant 



of the order of a wavelength of light. A local distortion of the periodicity in such structures 
can give rise to local photon states with frequencies within a gap between different photonic 
bands 0,0. 

Inspite all the obvious arguments regarding infeasibility to create local states of transverse 
electromagnetic field in convential crystals with microscopic defects, we show in this Letter 
that it is, nevertheless, possible. Our idea is to consider an interaction of an impurity with 
electromagnetic waves in a region of a polariton resonance, where light propagates as a 
mixed states with other excitations. Scattering of extended polaritons due to impurities was 



studied in Refs. ||13Hl5f . We consider local polariton states , which appear, due to defects, 
inside a forbidden gap between different polariton branches (reststrahlen region). We show 
that dipole active isolated defects can form coupled states with polaritons. If a frequency 
of this complex falls within the gap both electromagnetic and material components of these 
states are confined around the impurity, so that they can be called local polaritons. 

In this Letter we consider a simple model of a scalar harmonic excitation of a crystal 
(polarization wave) interacting with the electromagnetic field, which is also treated in the 
scalar approximation. A single defect embedded in the structure is assumed to affect the 
diagonal part of the corresponding Hamiltonian only, which means that the defect and host 
atoms differ only in their individual characteristics (e.g., masses for phonons or single atom 
excitation energies for excitons). Atom polarizability of the defect is assumed to be the same 
as for host atoms. This assumption makes our model of the defect considerably different 
from that considered in Ref. [|14[|, where defects were assumed to differ from host atoms only 



in their polarizabilities. We will show that effects due to elcetrical defectiveness are negligible 
comparing to those resulting from difference in mechanical (or electronic) characteristics. 

The general results presented here do not depend upon the specific nature of the crystal 
excitations (which we will call polarization waves) and when discussing them we will refer 
to both excitons and to phonons. We believe that this model reflects the most important 
characteristics of the phenomenon and will serve as a good starting point for more realsitic 
consideration. 

The electric component, E, of the transverse electromagnetic field interacting with scalar 
classical polarization waves obeys Maxwell's equation 
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c 2 dt 2 y ' c 2 dt 2 ' y ' 

where c is the velocity of an electromagnetic wave in the medium. The term on the right 
hand side of Eq. (JJ) is responsible for the interaction with lattice excitations due to their 
polarization, P. The equation of motion for the polarization waves can be written as 

d 2 P(r) d 2 P(r) d 2 

This equation describes one branch of excitations of a crystal with a single impurity at a site 
r interacting with the electric field E. The last term at the left side of Eq. represents 
a local perturbation due to the impurity, where the parameter 7 accounts for a difference 
between an impurity and host atoms. For phonons, 7 is the relative difference in masses, for 
excitons, it is the relative difference in atom's excitation energies. The dynamical matrix 
of a crystal, L(r — ri), which describes an interaction between atoms at different sites, is 
assumed to be unaffected by the impurity, parameter d is an effective coupling parameter of 
the excitations with the electromagnetic field. 

Polaritons arise as the simultaneous solutions of Eqs. (^) and (0). In order to derive an 
equation describing the electromagnetic component of polaritons, one needs to eliminate the 
polarization, P(r), from Eqs. ([]]) and (^). This can be done by using the Green's function 
of Eq. @, which has the form 



G(r, r') = G°(r - r') - 7 ^ 2 G°(r - r )G°(r - r^M -1 . (3) 

Here G°(r — r') is the Green's function of the polarization waves in a system without defects. 
Zeroes of D(u>), 

determine a frequency of local excitations in the absence of the interaction with the electro- 
magnetic field. In Eq. (|4]) summation over the first Brillouin zone is assumed, the function 
cjp(k) is an initial dispersion law for the polarization waves. In the long wavelength approx- 
imation this function is taken in the form 

u 2 p (k)=u; 2 -A 2 (ka) 2 , (5) 

which assumes a negative dispersion. The parameter A in Eq. (|5p gives the bandwidth of 
the spectrum of polarization waves. 

Using the Green's function (Rfi, one can obtain the equation for the spatial Fourier 
transform of electromagnetic field, E(k), in the form 

^ 2 E e kk>)£ kl = <?k 2 E k , (6) 

ki 

where the dielectric constant, ekki, is determined by the expression 

L d 2 \ _ rf^exp^k-kQro] 1 

£kkl I u? - u*(h) ) ^ 1 '[u, 2 - u$(k)][«? - <"2(*i)] D(uY [ ) 
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The diagonal part of the dielectric constant corresponds to the interaction of electromagnetic 
waves with the excitations of the ideal structure and in the absence of the defects it would 
lead to the usual polaritons. The second, non-diagonal, term in Eq. ([7|) reflects the effect 
of the impurity. Its additional pole due to localized excitations of the crystal will be shown 
to bring on local polaritons. Eq. (M) can be transformed into the following form 



where 



(u 2 - ujf)(u 2 - uj 2 p ) - d 2 u 2 \ E k = 7 dVexp (-ikr )r}{E)D{u), (8) 

^ E k exp (-ikr ) 



k ^ 2 - 

and ui(k) = ck is the dispersion law of electromagnetic waves in the crystal in absence of the 
interaction with polarization waves. This is a uniform integral equation with the parameter 
rj(E) depending upon all Fourier components, E k , of the electromagnetic field. In order to 
find the electromagnetic field from this equation one has to transform it into the form 

Parameter f](E) has a non-zero value only if the frequency satisfies the equation 



Solutions to Eq. ([TOf) determine new eigenfrequencies of the system which arise due to 
the impurity. We show that corresponding eigenmodes are localized in the vicinity of the 
impurity. The parameter f](E), though it is a function of E, is a free parameter determined 
by initial or boundary conditions for these frequencies. It has a simple relationship with the 
polarization on the impurity site: 

P(r ) = A7rd 2 7 1 (E)D- 1 (uj). (11) 

If the frequency does not satisfy Eq. (|10D , the parameter rj(E) is equal to zero and Eq. 
(§) describes ordinary polaritons, which are not affected by the defect. As it follows from 
Eq.flTTD the polarization of the defect in this case is also equal to zero, which means that 
the defect remains unexcited at all. 

Eq. ( |I0D gives rise to real eigenmodes if the expression in the sum in this equation does 
not have poles for real to and k. It occurs only if uj belongs to the gap between different 
polariton branches. This gap can exist if polarization waves have a negative dispersion at 
small k, which we assumed in Eq. (|5]). Then Eq. @ describes two branches of polaritons 
with a gap 5 P = u 2 — UJ\ between them. Here lo\ = u> 2 + d 2 , and uj\ ps uj 2 — Au d(a/c) is 
determined by the maximum of the lower polariton branch. This maximum is an important 
general feature of the polariton dispersion law in the presence of a bandgap. It has a 
topological origin and it always exists if polarization excitations have a negative dispersion 
in the long wavelength limit. We show that this is responsible for the significant difference 
between properties of states within the bandgap in our system and those within bandgaps 
caused by periodicity. 
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For frequencies far from the polariton region, Eq. (|T0| ) reduces to the well known form 
of Lifshitz's equation D~ 1 {oj) = 0, for the frequency of local phonon states. It is 

interesting to note, however, that for frequencies outside of the polariton bandgap, Eq. fllCf ) 
gains an imaginary part due to poles for real u, which describe a radiative decay of local 
vibrations or excitons. 

It is readily seen that for the frequencies inside the polariton gap the solution to the 
Eq. (Eft has the form E oc E (u)exp[— | r — r where the localization length, I, 

can be derived from the polariton dispersion law. Near the lower boundary of the gap, its 
dependence upon frequency has the form 



/ 1 oc \Juj 2 — Ui/aA, 

and near the upper boundary, 



7-1 ^2 / 2 9 



The frequency dependence of the localization length according to these expressions is asym- 
metric. It falls off rapidly when the frequency moves away from the lower boundary and 
increases slowly when the frequency approaches the upper boundary. The amplitude of the 
field E Q is given by the expression 



co 4 a 2 



E (u) » 7 P(r„) - (12) 

Ac 2 \ Ld 2 — (jjf 



where we have used Eq. ( |TTD in order to replace parameter rj(E) by -P(ro). 

The frequency uo and the parameter 7 in this expression are connected to each other 
by dispersion eq. (|TD|). This equation containes the sum over all wave vectors including 
those from the edge of the Brilluen band and from both polariton branches. At small k 
these branches are reduced into practically pure phonon and photon states. Photons with 
wavelengths comparable to a lattice parameter can not be treated in the framework of 
macroscopic theory used in this paper. However, it can be shown that contribution from 
photon-like branch in the region of large k is negligible (one can see this merely by neglecting 
the coupling parameter d in eq. (0)); and the only states contributing into the sum in eq. 
(P~0|) at large k are pure phonons. This elimination of photon contribution at the edge 
of Brilluen band is due to the mechanical nature of the defects under consideration. If 
one considers a defect atom with the polarizability different from that of host crystals the 
numerators in eq. (|10|) would have a form of u 2 — <Ji. Consequently, the major contribution 



to the sum at large k would arise from pure photon states. The sum with pure photon 
propagators, which arizes in this situation, is still assumed to be restricted by the first 
Brilluien band, since higher wave numbers can be thought as included into dynamical matrix 
L(r — ri). Local phonon and exciton states in one- and two-dimensional systems exist for 
an arbitrarily small value of 7. In 3-D systems, such states may appear only if 7 > 7 cr . 
This difference is due to the behavior of the corresponding integral for uj near the lower 
boundary of the gap: in 1-D and 2-D the integral diverges whereas in 3-D it has a finite 
value, which determines the threshold for 7. In the situation considered in this paper, the 
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sum in Eq. ( pTO|) diverges near the lower boundary of the gap u\ even in a 3-D situation. 
This singularity results from the maximum of the lowest polariton branch at to = oji. In 
order to show how this singularity occurs let us rewrite the sum in Eq. fllOD as follows 



El " ,7"" da) 
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[u 2 - nl(k)] [u 2 -n 2 2 {k) 



where Q\(k) and ^(k) are dispersion laws for lower and upper polariton branches respec- 
tively. Near the lower edge of the gap, fli {k) 2 can be presented as uj\ — b{k 2 — k 2 ^) 2 . Since the 
corresponding integral has a pole of second order at k = k m , it diverges when uo approches 

The sum in Eq. (|IUD cannot be evaluated in the long wavelength approximation since 
details of the spectrum of polarization waves for wave vectors near the boundary of the 
Brillouin zone contribute considerably to it. In the vicinity of the frequency u x , however, 
the only singular part of this sum becomes important. Since the singularity occurs at the 
wave number k m ~ (auod) / (cA) , which is greater than the resonance wave number (uq/c), 
but still far from the boundaries of the Brillouin zone, the singular part can be evaluated in 
the long wavelength approximation. 

Neglecting the final contributions to the integral, one can obtain the following relation 
between the frequency and the parameter 7 



^-"?~7^- (14) 

This equation shows that in the harmonic approximation, local polariton states arise for any 
value 0/7 in 3-D systems. This result demonstrates the difference between the two kind of 
gaps: those created by periodicity and by resonance interaction. This difference is due to 
that the density of states of excitations near the boundary of the first kind of gaps turns to 
zero, while in our case it diverges when the frequency approaches the lower boundary of the 
gap. 

According to Eq. ( fL4[) when 7 increases the frequency of the local state tends away from 
the boundary, inside the gap. This guarantees the existence of solutions at finite distance 
from the boundary at any finite value of 7 up to some 7 cr when the frequency becomes equal 
to the upper boundary of the gap oj%. For larger 7, the frequency runs out of the gap and real 
valued solutions of Eq. (|10D dissapear. The upper critical value, 7 cr , can be calculated from 
Eq. ([TOD with u = UJ2. The leading contribution to the integral at this frequency originates 
from k at the boundary of the Brillouin zone. In this situation the integral reduces to the 
corresponding integral in purely phonon or exciton problem. However, it is interesting to 
note that the shift of local frequencies from uj 2 into the gap caused by a small deviation of 
7 from 7 cr is strongly affected by the interaction with the electromagnetic field. The value 
of such a shift, 5u> 2 = uj\ — u 2 , can be described by the following expression 



(^7 2-K 2 0jla? rr- g 7l 2 UJ 
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cx 7 



where £7 = 7 cr — 7. The first term in this equation determines a; (7) near the upper boundary, 
which is similar to the result obtained for local mechanical vibrations ||: 002 — w oc 5~f 2 . 
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However, the coefficient of this term is extremely small so that, in the immediate vicinity 
of the upper boundary we have a crossover to the dependence u 2 — oj oc 5j, determined by 
the second term on the right hand side of Eq. (|1^). Hence, polariton effects result in a 
faster, linear shift of the local frequency from the boundary in comparison with quadratic 
"phonon-like" dependence. 

States arising in the immediate vicinity of a boundary of a gap are sensitive to different 
factors omitted in our model, such as anharmonicity or interactions with other subsystems. 
This is true for our model as well as for local excitations with regular behavior near the 
boundary of the gap. In general these factors result in the relaxation of all modes involved. 
This has a two-fold effect on the local states. First, they become quasistationary with a 
finite lifetime. This time, however, is usually several orders of magnitude greater then the 
oscillation period of local modes (e.g., for hydrogen- localized vibrational modes in GaAs this 
ratio is ~ I0 3 |16|]). Since the relaxation of local polaritons is mainly determined by the 
decay of phonon or exciton components, one can expect their lifetime to be approximately 
of the same order of magnitude. Second, due to relaxation processes, the boundary between 
the continuous spectrum and the gap acquires a finite width, T. This fact does not affect 
states deep inside of the gap, yet it influences states in the vicinity of the band edge. Due 
to the width T, these states appear on the background of the continuous spectrum and they 
cannot then be considered as truly local states. Therefore, only states which are father from 
the boundary than its relaxation width, T can be considered to be local. This increases the 
threshold for regular local states in periodic systems and sets up a finite threshold for local 
polaritons. This relaxation-induced threshold can be estimated from Eq. (|T4]) provided 
that the threshold is small enough that calculations of the integral in the main singular 
approximation still valid: 

„ a 2 rF~ 

(16) 

Numerical estimates show that for excitons with high frequencies and small damping it is 
possible to have for / ~f cr ~ I0~ 2 , where we have used parameters for exciton polaritons in 
KI from Ref. ||17|| . We conlcude, therefore, that one could observe local exciton-polariton 




states with frequencies within the polariton gap, induced by impurities with trap depth 
much smaller than for ordinary local excitons. 

Even if genuinely local states do not arise at small values of 7, this does not mean that the 
singularity in the Eq. (|10D has no effect. It brings on so called "virtual" or quasilocal states, 
which are also well known for phonon or exciton excitations. Even though these resonances 
are not real eigenstates of the system, they can considerably affect thermodynamical or 
optical properties of the material. In the situation considered, one can expect that these 
resonance states could influence the shape of the absorption edge in the vicinity of the lower 
boundary of the polariton gap. 

In this Letter we have shown that regular polar crystals with impurities can be used to 
create local states of polaritons with frequencies inside the polariton gap. Such local states 
in the case of a low concentration of impurities can manifest themselves experimentally in 
the frequency dependence of the transmission coefficient inside the forbidden gap. A similar 
effect has been observed for a single defect in a thin layer of a photonic crystal JTTJ|. One 



should expect a peak in the transmission when the frequency coincides with the frequency 
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of a local state. It occurs because an exponential tail of the incident field increases in 
the vicinity of the impurity if the frequency of the field coincides with the eigenfrequency 
of the local polaritons. These localized states provide one with as yet an unexplored yet 
opportunity to trap the transverse electromagnetic field and to produce narrow optical filters. 
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